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£N) ' Abstract 

| We consider the Vlasov-Einstein system in a spherically symmetric 

■ setting and prove the existence of static solutions which are asymptot- 

OO ' ically flat and have finite total mass and finite extension of the matter. 

Among these there are smooth, singularity-free solutions, which have 
a regular center and may have isotropic or anisotropic pressure, and 
f~>) ■ also solutions, which have a Schwarzschild-singularity at the center, 

(f") \ The paper is an extension of previous work, where only smooth, glob- 

Os • ally defined solutions with regular center and isotropic pressure were 

q ' considered, cf. || 

cr 

bjr)! 1 Introduction 

> ■ 

• rH , It is well known that the only static, spherically symmetric vacuum solutions 

r> | of Einstein's field equations are the Schwarzschild solutions which possess 

a spacetime singularity (or are identically flat). In the present note we 
couple Einstein's equations with the Vlasov or Liouville equation for a static, 
spherically symmetric distribution function / on phase space, describing an 
ensemble of identical particles such as stars in a galaxy, galaxies in a galaxy 
cluster etc.. This results in the following system of equations: 



^— ■d x f-yjl + v 2 f i'~-d v f = 0, x,v£M 3 , r:=\x\, 



e" 2A (2rA'-l) + l = %i:r 2 p, 
e- 2A (2?y + l)-l = 8vrr 2 p, 
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where 



p(x) = p(r) 
p{x) =p(r) 

Here the prime denotes derivative with respect to r, and spherical symmetry 
of / means that f(Ax,Av) = f(x,v) for every orthogonal matrix A and x,v € 
K 3 . If we let x = r(sm9cos4>,sm6sm(f>,cosd) then the spacetime metric is 
given by 

ds 2 = -e 2 ^dt 2 + e 2X dr 2 + r 2 {d6 2 + sin 2 ## 2 ). 
Asymptotic flatness is guaranteed by the boundary condition 

lim A(r) = lim u(r) = 0, 

r— >oo r-^oo 

a regular center at r = is guaranteed by the boundary condition 

A(0) = 0. 

We refer to || on the particular choice of coordinates on phase space leading 
to the above formulation of the system. It should be pointed out that 
the above equation for / is only equivalent to the Vlasov equation if / is 
spherically symmetric. 

A brief overview of the literature related to the present investigation 
seems in order. In the initial value problem for the corresponding time 
dependent system is investigated and the existence of singularity-free solu- 
tions is established for small initial data. In it is shown that the classical 
Vlasov-Poisson system is the Newtonian limit of the Vlasov-Einstein system. 
The existence of globally defined, smooth, static solutions with a regular cen- 
ter, finite radius, and finite total mass is investigated in Q; that investigation 
is restricted to solutions with isotropic pressure. For the Vlasov-Poisson sys- 
tem the existence of spherically symmetric steady states is established in || , 
and we briefly recall the general approach followed there: Taking the distri- 
bution function / as a function of the local energy and the modulus of the 
angular momentum, which are conserved along characteristics, the problem 
is reduced to solving a nonlinear Poisson problem. For certain cases corre- 
sponding to the so-called polytropes it can then be shown that this leads to 
solutions with finite radius and finite total mass, cf. also [||. Cylindrically 
symmetric, static solutions of the Vlasov-Poisson system are investigated in 
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+ v 2 f(x,v)dv, 

x-v\ 2 dv 
f{x,v) 



R 3 V r 



Vl + v 2 ' 
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[0. A system in some sense intermediate between the Newtonian and the 
general-relativistic setting is the relativistic Vlasov-Poisson system, where 
similar results as for the nonrelativistic version hold, cf. JjJ]. In (II]] it is 
shown that coupling Einstein's equations with a Yang-Mills field can lead 
to static, singularity-free solutions with finite total mass. Static solutions of 
general-relativistic elasticity have been studied in |J. For further references, 
especially on the astrophysical literature, we refer to fTpf . 

The present investigation proceeds as follows: In the next section we 
give two conserved quantities, which correspond to the local energy and 
angular momentum in the Newtonian limit, and reduce the existence prob- 
lem for static solutions of the Vlasov-Einstein system to solving a nonlinear 
system of ordinary differential equations for A and fi. The fact that the 
distribution function / indeed has to be a function of the local energy and 
angular momentum — usually referred to as Jeans' Theorem — is rigorously 
established for the Vlasov-Poisson system in Q but not known in the present 
situation. In other words, it is not clear whether our ansatz for / covers all 
possible spherically symmetric, static solutions. The existence of solutions to 
the remaining system for A and ji is investigated in the third section. There 
we dispose of the restriction that / depends on the local energy only, which 
was used in || and lead to solutions with necessarily isotropic pressure. 
Nevertheless, the arguments used in the present, more general situation are 
considerably simpler than the ones used in ||; a key element in the proof 
that the solutions exist globally in r is the Tolman-Oppenheimer-Volkov 
equation. In the fourth section we show that — for an appropriate ansatz 
for /, corresponding to the polytropes in the classical case — the solutions 
obtained by our approach have regular center, finite mass, and finite radius. 
This is done by treating the Vlasov-Einstein system as a perturbation of 
the Vlasov-Poisson system and using the criteria for a finite radius from the 
classical Lane- Emden- Fowler equation. In the last section we prove the exis- 
tence of solutions with a Schwarzschild-singularity at the center, surrounded 
by Vlasov-matter with finite radius and finite total mass. Surprisingly, the 
latter features of these solutions can be obtained easily and without the 
perturbation argument used in the nonsingular case mentioned above. 
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2 Conserved quantities and reduction of the prob- 
lem 



As explained in the introduction, a key point in our investigation is to reduce 
the full system to a nonlinear system of ordinary differential equations for A 
and p. This is achieved by the ansatz that the distribution function depends 
only on certain integrals of the characteristic system. In the coordinates used 
above, the characteristic system corresponding to the above equation for / 
reads 

v 

-\Jl + v 2 p'(r) -; 

r 

note that these equations are not the geodesic equations. One immediately 
checks that the quantities 

E:=e^y/l + v 2 , F:=x 2 v 2 - {x-vf = \x x v\ 2 

are conserved along characteristics; in Sect. 4 the relation of E to the classical 
local energy will become apparent, F can be interpreted as the modulus of 
the angular momentum. If we take / to be of the form 

f(x,v) = *(E,F) 

for some function the Vlasov equation is automatically satisfied. Inserting 
this into the definitions of p and p we obtain the following nonlinear system 
for A and p: 

e" 2A (2rA' - 1) + 1 = 8irr 2 G$(r,p), (2.1) 
e~ 2A (2ry + 1)-1 = 8wr 2 H^(r,p). (2.2) 

Here 

27r roo r r 2 (e 2 -l) ,2 

G*(r,«) = — / / $(e^ r h,F) -== = dFde, 

r l J i Jo y/e 2 -l- F/r 2 

H^(r,u) = -^J^ <S>{e^e,F)^e 2 -\-F/r 2 dFde 

for r, u > which follows by a simple transformation of variables in the in- 
tegrals defining p and p. These latter quantities are then given by 

p(r) = Ga>(r,n(r)), p(r) = H<s>(r, p(r)) , 



x = 

V = 
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once a solution of the system (|2.1|), is known to exist. 

Of course one has to make some assumptions on $ in order to obtain 
existence results for the above system and to investigate the properties of its 
solutions. The assumption used below is certainly not the weakest possible, 
but on the one hand it is sufficiently general to encompass a large class of 
examples, leading to quite different classes of solutions, and on the other 
hand the main ideas of the proofs are not buried under technicalities. We 
require $ to satisfy the following 

Assumption on $: 

$(E,F) = <f>(E){F-F ) l + , E>0, F>0, 

where Fq>0, l>— 1/2, and € L°°(]0,oo[) is nonnegative with <p(E) = 
0, E > E , for some E > 0. 

Under this assumption the functions G$ and H$ are easily seen to take 
the form 



G^r,u) = Cl r 2l e-^ 




H*(r,u) = ^L_r 2l e-^ u h^ (e u ^l + F /r^ , (2.4) 

where 

roc , 

g^t):=J <P(E)E 2 (E 2 -t 2 ) l+ 2dE, (2.5) 
fy(t):=y 0(E)(E 2 -t 2 ) l+ idE (2.6) 

for t > 0, and 

ci :=2tt / ds. 

Jo Vl — s 

3 Existence of solutions 

As a first step we show that the functions g$ and defined by fl2,5] ) and 
( |2.6D are well behaved for functions <f> as considered in the assumption above. 
More precisely: 
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Lemma 3.1 Eqns. (2.5) and $2.(\ ) define decreasing functions g^,,h^€z 
C 1 (IR + ) and 

tij(t) = -(2l + 3)t J (f}(E)(E 2 -t 2 ) l+1 / 2 dE, t>0. 

Proof: Obviously, the integrals defining and exist. Lebesgue's 
dominated convergence theorem implies that the functions g^ and are 
continuous. For t > and At > such that t — At > we have 

±(h^t-At)-h^t)) = ^J t [^(E)(E 2 -(t-At) 2 ) l+V2 dE 

+ / ^)^({E 2 -{t-At) 2 ) l+z l 2 -{E 2 -t 2 ) l+ ^ 2 )dE 
=:h + h- 

Obviously, I\ —* as At — > 0. The term /2 has a limit as At —* by Lebesgue's 
theorem. Thus, the function /i^ is left-differentiable with 

- h^t) = -{2M)tJ <p(E)(E 2 -t 2 ) l+1 / 2 dE. 

Again by Lebesgue's theorem this function is continuous, and a continuous 
and continuously left-differentiable function is continuously differentiable. 
Similarly, we see that 

- g4t) = -(2l + l)tj t <P(E)E 2 {E 2 -t 2 ) l - 1 l 2 dE, 

which is again continuous. The derivatives are negative, and thus g^ and hj, 
decrease. □ 

The regularity of the functions g^ and being established, we can now 
prove local existence and uniqueness of A and [i. Note that in the last section 
we will need to pose initial data at some ro > which is why we include this 
case in the following results. 

Theorem 3.2 Let $ satisfy the general assumption and let G$ and H§ 



be defined by Eqns. $U\) , $A), ( Hjy - Then for every r > and 

Ao,/xoSM with Ao = z/ro = there exists a unique maximal solution A,//€ 
C x ([ro,R\), R>r , of the Eqns. $L%) and ( ggj with 

A(r- ) = A , (j,(r ) = no- 
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Proof: Due to the regularity of the right hand sides in Q2.1|), ( |2.2| ) — cf. 
Lemma |3.1| — the result is of course trivial in the case tq > 0. In the case 



ro = one has to deal with the singularity of the Eqns. ( |2.1| ) and (2.2) at 
r = 0. Using the boundary condition at r = 0, Eqn. (|2.1| ) can be integrated 
to give 

e -2A(r) = 1 _^ f r s 2 G^(s,fi(s))ds. (3.1) 
r Jo 

If we insert this into Eqn. (^^) we obtain an equation for fi alone: 

(3.2) 

Integrating Eqn. ( |3.2[ ) subject to the initial condition fj,(0) = /jlq we obtain 
the following fixed point problem for \x: 

H(r) = (TfJt)(r), r>0 

where 

(Tfj)(r):=W) + 



4tt 



sH$(s,fj,(s)) + — J a 2 Gq,(a,fi(a))da^j ds. 
A lengthy but straight forward argument shows that T maps the set 

M:={/i:[0,5] — H| fi{0) = no, W<mW<W + 1, 

— [ r S 2 G*(s,^s))ds<l, rG[0,5]\ 
r Jo 2 > 

into itself and acts as a contraction with respect to the norm ||-||oo, pro- 
vided 5 > is chosen small enough. This gives the existence of a solution of 
Eqn. ( |3.2| ) on the interval [0,5]. On this interval we define A by Eqn. ( |3.1| ) 
and obtain a local solution A,/i G C 1 ([0,<5]) of fl2,l| ), (|2.2|) . Obviously, the 
boundary conditions at r = are satisfied, and the solution is unique. □ 

In order to show that the above solutions actually extend to r = oo we 
need a relation which is known as the Tolman-Oppenheimer-Volkov equation 
in the context of the general relativistic description of spherically symmetric 
fluid balls, cf. pi, Eqn. 6.2.19]. 
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Lemma 3.3 Define 



PT{r):= \ [ 
2 Jn 





XXV 


fn 3 


r 



f(x,v) 



dv 



and let X,fi be a solution of the system \2.\ ), j2.2\j on the interval [ro,i?[. 
Then px has the form 

p T (r) = (l + l)p(r) + ^F r 2l - 2 e-( 2l+2 ^ kj, [e^ y/l + F Q /r^ , 

where 

jfe^(t):=y (j){E){E 2 -t 2 ) l+l,2 dE, t>0, 
defines a decreasing C 1 -function, and 



p'(r) 



-fjf(r) {p{r) + p{r)) - - (p(r) -pr{r)) , r e]r ,R[. 
r 



Proof: The formula for pt is obtained by a transformation of variables, 
and the regularity and monotonicity of k^ follow as in Lemma 3.1. Using 
(|2.4j) , (2.6), and Lemma one obtains the relation 

21 

p'(r) = — p(r) — (21 +4)/i (r)p(r) 
r 

-q (l + F /r 2 ) r ^ e -(2i+2)^(r) ^ ^(r)y i + F()/r 2j ^ (r ) 



+c l r 2l -*F Q e- { - 2l+2 ^k <t) e^Jl + F /r 2 



which, after some further calculations, leads to the Tolman-Oppenheimer- 
Volkov equation. □ 

Remark: The quantity pt is the tangential pressure generated by the phase 
space density /, as opposed to the radial pressure p. Note that if Fq = then 
p and pt differ only by a multiplicative constant, and if Fq = 1 = 0, i. e. the 
phase space distribution is independent of the angular momentum, then 
P = Pt, i- e. the solution has isotropic pressure. 

Theorem 3.4 Let <3? satisfy the general assumption and let G$ and H§ be 



defined by Eqns. ftjUj ), $Q0, < $A )- Then f or ever V r > and A > 0, 

HoGlR with Aq = i/ro = there exists a unique solution X,^G C 1 ([ro,oo[) 



of the Eqns. ft2J\ ) and $2J\ ) with 

A(r ) = A , /i(r ) = Mo- 
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Proof: Let A,^G C 1 ([ro,-R[) be the maximal solution to (|2.1|), ( |2.2[ ) which 
exists by Thm. |3.2| . We have the following equations on [ro,i?[: 

e~ 2A (2rA / -l) + l = 87rr 2 p(r), (3.3) 

e" 2A (2r^' + l)-l = 87rr 2 p(r), (3.4) 
2 

p (r) = (r) (p(r)+p(r)) - - (p(r) — pr(»")) • ( 3 - 5 ) 

Note that — as opposed to the case considered in || — the functions p and 
p are not necessarily strictly decreasing on their support so that we cannot 
write p as a function of p, and we cannot use the analysis in Q directly. 
Nevertheless, the proof given here is somewhat simpler than the one given 
in H . If we integrate Eqn. ( |3.3D we obtain 



e -2A(r) = 1 _8vr r 2 p(s)ds _I r , o(1 _ e -2A 0)> (3 g) 

' Iro r 



r 



Inserting this into Eqn. (3^) yields 

p'(r)=47re 2X{r) (p(r)+w(r)) (3.7) 

where 

. W; =-L(£ S VW* + ^^). (3.8) 

Finally, by adding Eqns. ( |3.3| ) and ([TJ|) we have 

(A' + p') (r) = 47rre 2A(r) (p{r) +p(r)). (3.9) 

We now wish to establish a differential inequality for the quantity e x+tJ, (p + 

w), which will allow us to bound A and p and conclude that R = oo. Using 
Eqns. dp]), Q, and Q we see that 

2p 2pr 3w | p 



(e x+ »(p+w)y (r) = e x+fl f-y 



\ r r 

note that the terms which were dropped are indeed negative by the as- 
sumption Ao > 0. Now assume that R<oo. Without loss of generality, we 
may assume that there exists r\ €]ro-R[ and C>0 such that w(r)>C for 
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r£ [ri,i?[, since otherwise the solution is trivial. The fact that p is increas- 
ing and the functions and are decreasing implies that p and pt and 
thus also p/r and pt/t are bounded on the interval [ri,i?[. Thus, we can 
continue the above estimate to obtain 

(e x+ ^ (p + w)) ' (r) < Ce x+tI < Ce x+fl (p(r) +w(r)). 

This implies that 

e A(r)+„(r)< < c re[n,i?[. 

Since A(r) > by ( ^ ) and /x(r) > /Xq by monotonicity, this implies that A 
and p are bounded on [?"o,i?[ which is a contradiction to the maximality of 
R. Thus i? = oo, and the proof is complete. □ 



Remarks: 

1. If r >0 then X,p G C 2 ([r , oo[) and C 1 ([r ,cx)[). 

2. If r = F = then A,/^G C 1 ([0,oo[) n(7 2 (]0,cx)[) with A'(0) = //(0) = 
and C 1 (]0,oo[). If in addition Z = or Z > 1/2 then A,/z<E 
C 2 ([0,oo[)nC 2 (IR 3 ) and C^^oo^nC 1 ^ 3 ) with p'(0) =p'(0) = 
where functions in r are identified with the corresponding radially 
symmetric functions on IR 3 . Thus we see that in the case r$ = the 
solutions are as regular at the center as anywhere else if various pa- 
rameters are chosen appropriately. 

3. The phase space density / := $>(E,F) is a solution of the Vlasov equa- 
tion in the sense that it is constant along characteristics. If $ is con- 
tinuous or continuously differentiable then in addition / has the same 
regularity, and in the latter case satisfies Vlasov's equation classically. 



4 Singularity-free solutions with finite mass and 
finite radius 

In this section we are interested in smooth, singularity-free solutions with 
with a regular center. Thus, we consider the boundary condition A(0) = 0. 
Also, we set Fq = in this section, the case -Fo>0 will play its role in the 
next section. 
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In order to decide whether a solution obtained in the previous section 
has finite total mass or whether p(r) vanishes for r large, one has to have 
rather detailed information on the behaviour of the function p. Due to 



the complexity of Eqn. (3.2) we have not been able to decide these ques- 
tions directly, even for simple examples of $ and without dependence on F. 
However, it is possible to show that the solutions obtained above converge 
to solutions of the corresponding Newtonian problem as the speed of light 
tends to infinity. It is then possible to use the results on finite mass and 
finite radius which are known in the Newtonian case for the so-called poly- 
tropes to obtain solutions with the same properties for the Vlasov-Einstein 
system. 

To carry out this program we introduce the parameter 7 := \ where 
c denotes the speed of light, define v := ^p, and recall from 0] that the 
Vlasov-Einstein system with 7 inserted in the appropriate places reads 



■d x f- x /l + 7 v 2 v'--d v f = 0, 
r 



where 



e - 2A (2rA'-l) + l = 8irjr 2 p, 
e _2A (2rz/ + 1/7) - I/7 = 8vr7r 2 p, 



p(x) := / J 1 + jv 2 f (x ,v) dv , 

dv 



The conserved quantity E now becomes 



and F remains unchanged. In order to obtain the correct limit as 7^0, 
we have to rewrite our ansatz for the distribution function in the following 
form: 

f{x,v) = 4> ^1 + jv 2 e iv( - r) - F l . (4.1) 

As above the Vlasov-Einstein system can then be reduced to a single equa- 
tion for v alone, namely 

i/(r)= (4.2) 
1 8 , r rr^T (ir 2l+1 h^(v(r)) + ± f ' s 2l+2 g^(u(s))ds) , 
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where 

roo / i i \ £+1/2 

^ 7 (n):=qe- 2 (^y^^)(l+ 7 ii;) 2 U(l + 7 i?) 2 --e 2 ^J dE, 

7 7 7 

(4.3) 



(4.4) 



so that 

„2Z„ „/„A „2J 



p(r)=r%, 7 (i/(r)), p(r) = r 2i /i 0i7 (z^(r)). (4.5) 
For the Newtonian case the corresponding ansatz 



,2 



f( x , v ) = <j)\- + U{r)\F l 



reduces the Vlasov-Poisson system to the equation 



U'{r) = % I s 2l+2 g (U(s))d S , (4.6) 



r 2 Jo 



where 

go(u) := ci I " 0(B) (2(B - u)) l+1/2 dE. 



Assume that 0£L°°(IR) and Bo = 0, and fix uqKO. Clearly, the results of 
the previous section apply so that for every 7>0 Eqn. (4.2) has a unique, 
nontrivial, global solution with v(0)=vq. Let U G C 1 ([0,oo[) be the global 
solution of ( |4.6| ) with U{fS) = VQ. We shall prove that v converges to U as 
7^0, more precisely: 

Lemma 4.1 For every R>0 there exist constants C>0 and 7 o>0 such 
that for every 7E]0,7o] the solution v of \4-<\ ) with v(§) = vq satisfies the 
estimate 

\v(r)-U(r)\<C"f min{l+1/2 ' 1} , re [0,R], 
where U is the solution of ( \j.o\ ) with £7(0) = ^o- 

Proof: We have the estimate 

W'(r)-U'(r)\<h + I 2 + I 3 +h, 
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where 



47T 



J-l ■- 

h ■-- 
h--~ 



l-^lf s^g^{v{s))ds 



7r 2 ' +1 /^>(r)), 



1 



47T 

r 2 Jo 



rr 

Jo 



l-^lf s 2 ^ 2 g^{u{s))d, 

47T f r 

— s 2l+2 \ g<j)n {v{s))-go{y{s))\ds, 

i J 

4-tt r r 

J 4 == -^2 I s 2l+2 \go(v(s))-go(U(s))\ds. 

Since for r > we have v{r) > vq and since 4> is bounded with 4>(E) = for 
E 1 > 0, it is easily seen that 

g^{y{r))<C, h^{u{r))<C, r>0, 7 e]0,l], 

where the constant C depends only on <p and i>o- Thus, for R> fixed and 
r£ [0,R] we obtain the estimate 



1 



l-CR 2l + 2 ~f 



Ci? 2+1 7 < C 7 



where 7G]0,7o], and 7oG]0,l] is chosen such that 1 — CR 2l+2 ^o > 0- Simi- 
larly, 

/ 2 <c 7 

for 7 g]0,7o]. Next we estimate the difference ^^(^(r)) — go(v(r))\. We can 
restrict ourselves to the case v{r) < since otherwise this difference is zero. 
Thus, 

->Kr), 

7 

which implies that we can split the difference in the following way: 
\H</>( v ( r )) -9o{v{r))\ < C(Ji + J 2 + J 3 ), 

where 

1+1/2 

Ji := 



2(I+2) 7 ,_ 1 | r_^ (£ ; )(1 + 7£0 2/l (1 + 7£0 2_l^ ^ 

J — - — \7 7 / 

J 2 := / 7 <p(E)(2(E-u)) l+1 ^ 2 dE, 

J V 
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fOO 

el 



i 



i+7^r ( -(1+7^0 



2 _ i p 2 7 v 



1+1/2 



(2(E-u)) l+1/2 



dE. 



Now 



2(l+2)jv(r) _ i 

Ji < c 7 < c 7 , 



J 2 < c 



7 
7 



and 



(l +7 £) 2 -(1 + 7^ 



7 



J+l/2 



■ i/(r) <c 7 , 



{2{E-u)) l+1/2 



<C 
<C 



(1 + jE) 2 -1 



(1 + 7 E) 2 -1 



+ 



+ C 



1 1 \ '+1/2 

L {1 + lE f_i e 2ju 

7 7 



(2(^-1/)) 



1+1/2 



<C 7 + C 

<c 7 +c 7 



l + 2 7 z^(r) — e 



27^(r) 



7 

min{Z+l/2,l} 



+ 7^ 



7 

min{2+l/2,l} 



(2(^-1/)) 



2+1/2 



Finally, using the method of Lemma 3J. it is easily seen that goeC^M) 
and |<7o(w)| f° r u > ^o- Thus 



r C 
r z Jo 



21+2 



u{s)-U(s)\ds. 



If we combine all the above estimates we see that 



W{r)-U'{r)\<C^ 



min{Z+l/2,l} 



+ 



c_ 

r 2 Jo 



^21+2 



u(s)-U(s)\ds 



for rG [0,R] and 7 G]0, 7 o]. The assertion of the lemma now follows by the 
usual Gronwall argument. □ 

In order to exploit the above result, we have to restrict our investigation 
to such cases where finiteness of the radius and total mass is known in the 
Newtonian situation. Thus we consider 



^(E,F):=(-E)lF l ,E£M, F>0, 



(4.7) 
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where £;>0 and l> — 1/2 are such that k<3l + 7/2. Then we know from 
H 5.4] that the Newtonian potential U has a zero, which is also the radius 
where the density vanishes. Since by the assumption vq < the solution 
is nontrivial, U has to strictly increase so that there exists R > such that 
U (R) > 0. The above lemma then tells us that for all 7 > sufficiently small, 



v(R)>0 which by the definition of <3? and Eqn. (4.5) implies that p(r 



p(r) = for r>R. Thus we have proved the following theorem: 

Theorem 4.2 Let <E> be defined by Eqn. ( \4- 1 ) and take i/q < 0. Then for all 
7 > sufficiently small the corresponding solution v has a zero, the density p 



and radial pressure p as defined by Eqns. (4-5), (4-3 )> an d (4-4) have finite 
support, and 

poo 

0<4ir s 2 p(s)ds< 00 
Jo 

i.e. the solution is nontrivial with finite total mass. 



Remarks: 

1. In this section we have obtained static solutions of the spherically 
symmetric Vlasov-Einstein system with finite radius and finite mass, 
provided 7 > is sufficiently small. However, if /, A, v is such a solution 
then 7" 3/2 /(7" 1/2 -,7~ 1/2 -),A(7 _1 / 2 -),7^(7" 1/2 -) is a solution of the 
system with 7 = 1, which again has finite radius and finite total mass. 

2. If $ is defined by Eqn. ( p~7| ) then at least / e C 1 (supp/) , and fe 
COR 6 ) or /eC^H 6 ) if k>0,l>0 or fc>l,J>l respectively. The 
Vlasov equation then holds classically on supp/ or on IR 6 respectively. 

3. The finiteness of the total mass implies that A(r)^0 as r— >0 and v 
and p have a finite limit as r^oo. This means that the spacetime 
is asymptotically flat; the fact that the limit of p is not zero only 
corresponds to a rescaling of the time coordinate t. 

4. Our solutions are not only global in the coordinates which we used, 
but are singularity-free in the sense that the corresponding spacetime 
is timelike and null geodesically complete. 
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5 Solutions with a Schwarzschild singularity at the 
center 



In this last section we take $ as in Sect. 1, with Fq > 0, Eq = 1, and <j>(E) > 
for E < 1. This implies that 



G$(r,u) = H<s,(r,u) = if e u Jl + F /r 2 > 1, 



G$(r,u),#$(r,it)>0 if eVl + F /r 2 <l 



Consider the field equations ( |2.1| ) , ( |2,2[ ) . It is well kwown that in the vacuum 
case, i. e. if the right hand sides are identically zero, the asymptotically flat 
solutions are given by 

e Mr) = l _ 2 _M» j e 2X(r) = ( 1 _^.Y\ r>2M0} 

r \ r ) 

the so-called Schwarzschild metric; Mq > 0. This metric is investigated in 
probably every textbook on general relativity, it represents a spacetime sin- 
gularity which is hidden inside of an event horizon, a black hole. In passing 
we mention that the apparent singularity at r = 2Mq is only a coordinate 
singularity, i. e. in other coordinates the spacetime can be extended beyond 



this radius, which is also called Schwarzschild radius, cf. [12, Ch. 6]. In the 
present section we wish to construct static solutions of the spherically sym- 
metric Vlasov-Einstein system which have such a black hole at the center. 
Let us first see if or where the Schwarzschild metric solves our system 



(|2.l[) , d2.2|) . This is obviously the case if 



i-^/i+5>i, 



/- 



i. e. for r£ [r_,r + ], where 



Fo±JF$-WMgF 
r+ := . 

If we take Fq > 16Mq > we obtain 2Mq < r_ < r+, and we may set 

e 2M(D = 1 _^0 e 2X(r )= ( l _2AV- 
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and 

p(r)=p(r)=f(x,v)=0 
for 2Mq < r < r + . At r$ := r + we now pose initial conditions 

M(»"o) = Mo:=Jl -,A(r ) = A :=Wl , 

V r o V r o 



and solve ( |2.1[) , (|2.2| ) to the right of ro according to Thm. 3A; note that 
Ao > 0. In this way we obviously obtain a static solution of the spherically 
symmetric Vlasov-Einstein system, which coincides with the Schwarzschild 
solution for r<ro. The phase space density / is no longer given by one 
function of E and F because the vacuum solution is not consistent with 
our ansatz f = <&(E,F) in the region ]2Mo,r_[, but / still solves the Vlasov 
equation everywhere in the region ]2Mo,oo[ since no characteristics which 
carry mass can cross the region ]r_,ro[. 

What remains to be seen is that our solution is not vacuum everywhere 
for r > 2Mq, i. e. we really get Vlasov matter for r > ro, and that this matter 
outside ro again has finite radius and finite total mass. First note that 



<0 



which implies that p(r) > and p{r) > in a right neighborhood of ro- To 
show that for large enough values of r the density vanishes again, it is 
obviously sufficient to show that lim r ^ 00 ^(r) >0. To see the latter we first 
recall that by integrating Eqn. ( |2.1| ) we obtain 




e-W = 1 _MlO_I ro(1 _ e - 2 A 0) = x _ 2MW r > 2Moj 



r 

where 



rr 

M(r):=M + m(r):=M + Air / s 2 p(s) ds. 

Jro 

Inserting this into Eqn. (|2.2| ) we obtain 



r 

> M | m(r) 
- r (r-2M ) r 2 
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and integrating this yields 



rr d s f r m(s) 

fi(r)>no + M ^TT\ + ~^ ds 

Jr s(s — 2M ) J ro S z 



= w+1 ^ 1 _^£_ ln ^ 1 _^ ++ £=i(£) ds 



ro 



S 



Therefore 

. , f°°m(s) , 

'■ Jro S 2 



lim 



which proves that the density vanishes for r large enough. Note also that 
the inner boundary of the Vlasov matter satisfies the estimate 

ro> 4§o >4M °- 
We collect these results in the following theorem: 

Theorem 5.1 Let <I> satisfy the general assumption in Sect. 1 with Eq = 1 
and 4>(E) > for E>1, and let F > 16Mq > 0. Then there exists a static 
solution (/,A,/x) of the spherically symmetric Vlasov- Einstein system such 
that 

e 2 Mr) = 1 _2Mo ) e 2A W = / 1 _2Mor 1 ) 2Mo<r<r0) 
r \ r ) 



and 



p(r) =p(r) = f(x,v) = 0, 2Mo<r<ro orr>R, 
f R 

0<4tt p(r)dr <oo, 

Jro 



where 

F + v /f 2 -16M 2 F 
r := L_ >4 Mo 

and R>r . Furthermore, A,/xG C 2 (]2M ,oo[), C 1 (]2M ,oo[), and the 
spacetime in asymptotically flat. 

Remarks: 
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1. As pointed out above, the phase space density / is not given as a 
function of E and F globally for r > 2Mq. 

2. If we start at the center with a smooth static solution as obtained 
in Sect. 4 instead of the Schwarzschild singularity we can use the 
method of the present section to obtain static solutions of the fol- 
lowing kind: For < r < vq the nontrivial matter distribution is given 
by f(x,v) = &o(E), for ro < r < n we have vacuum, for r\<r<V2 the 
matter distribution is again nonzero and given by f(x,v) = $>i(E,F), 
and this procedure can be continued. The resulting solution is smooth, 
geodesically complete, asymptotically flat, with finite mass and finite 
radius, and consists of rings of Vlasov matter separated by vacuum. 
Contrary to our expectation it is not clear that such a solution violates 
Jeans' Theorem, which in spite of 1) might still be valid for smooth 
solutions. Note however that as the radius increases, the pressure can 
change from isotropic to unisotropic. 

3. In the Newtonian case one can construct analogous solutions, where 
the role of the Schwarzschild singularity is played by a point mass 
Mo situated at the center r = 0. It then turns out that the resulting 
solution does not violate Jeans' Theorem — in the above notation, r_ = 
in the Newtonian case. Thus one can at least say that the range 
of validity of Jeans' Theorem is larger in the classical case than in 
the general relativistic. One can also construct solutions of the type 
described in 2) for the Newtonian case. 

References 

[1] Batt, J.: Steady state solutions of the relativistic Vlasov- Poisson sys- 
tem. In Ruffini, R. (ed.): Proceedings of the Fifth Marcel Grossmann 
Meeting on General Relativity. 1988 

[2] Batt, J., Berestycki, H., Degond, P., Perthame, B.: Some families of 
solutions of the Vlasov-Poisson system. Arch. Rational Mech. Anal. 
104, 79-103 (1988) 

[3] Batt, J., Faltenbacher, W., Horst, E.: Stationary spherically symmetric 
models in stellar dynamics. Arch. Rational Mech. Anal. 93, 159-183 
(1986) 



19 



[4] Batt, J., Pfaffelmoser, K.: On the radius continuity of the models of 
polytropic gas spheres which correspond to the positive solutions of the 
generalized Emden-Fowler equation. Math. Meth. in the Appl. Sci. 10, 
499-516 (1988) 

[5] Kijowski, J., Magli, G.: A generalization of the relativistic equilibrium 
equation for a nonrotating star. General Relativity and Gravitation 24, 
139-158 (1992) 

[6] Rein, G., Rendall, A.: Global existence of solutions of the spherically 
symmetric Vlasov-Einstein system with small initial data. Commun. 
Math. Phys. 150, 561-583 (1992) 

[7] Rein, G., Rendall, A.: The Newtonian limit of the spherically sym- 
metric Vlasov-Einstein system. Commun. Math. Phys. 150, 585-591 
(1992) 

[8] Rein, G., Rendall, A.: Smooth static solutions of the spherically sym- 
metric Vlasov-Einstein system. To appear in Ann. de l'lnst. H. Poincare, 
Physique Theorique 

[9] Rendall, A., Schmidt, B.: Existence and properties of spherically sym- 
metric static fluid bodies with a given equation of state. Class. Quantum 
Grav. 8, 985-1000 (1991) 

[10] Shapiro, S. L., Teukolsky, S. A.: Relativistic stellar dynamics on the 
computer, I. Motivation and numerical method. Astrophysical Journal 
298, 34-57 (1985) 

[11] Smoller, J., Wasserman, A., Yau, S.-T., McLeod, J.: Smooth static 
solutions of the Einstein/ Yang-Mills equations. Commun. Math. Phys. 
143, 115-147 (1991) 

[12] Wald, R. M.: General Relativity. The University of Chicago Press, 1984 



20 



